Plaquette Boson-Fermion Model of Cuprates 
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The strongly interacting Hubbard model on the square lattice is reduced to the low energy Plaque- 
tte Boson Fermion Model (PBFM). The four bosons (an antiferromagnon triplet and a d-wave hole 
pair), and the fermions are defined by the lowest plaquette eigenstates. We apply the Contractor 
Renormalization method of Morningstar and Weinstein to compute the boson effective interactions. 
The range-3 truncation error is found to be very small, signaling short hole-pair and magnon coher- 
ence lengths. The pair-hopping and magnon interactions are comparable, which explains the rapid 
destruction of antiferromagnetic order with emergence of superconductivity, and validates a key 
assumption of the projected SO (5) theory. A vacuum crossing at larger doping marks a transition 
into the overdoped regime. With hole fermions occupying small Fermi pockets and Andreev coupled 
to hole pair bosons, the PBFM yields several testable predictions for photoemmission, tunneling 
asymmetry and entropy measurements. 



I. INTRODUCTION 

In 1987, shortly after the discovery of high temperature 
superconductivity in cuprates, AndersonEI proposed that 
the key to this perplexing phenomenon hides in the large 
positive Hubbard interactions in the copper oxide planes. 
Indeed, at zero hole doping, the Hubbard model captures 
the Mott insulator physics of the parent compounds e.g. 
La2CuC>4. The doped Hubbard model however has so far 
resisted a definitive solution, primarily because its spins 
and holes are highly entangled with no obvious small pa- 
rameter to separate them. Whether the Hubbard model 
even supports superconductivity without additional in- 
teractions remains a subject of controversy. Different 
mean field theories suggest conflicting ground state or- 
der parameters and correlations. Numerical methods-are 
restricted to finite clusters where hole pairing is founcoQ, 
but off-diagonal long range order has not been ascer- 
tained. 

This paper charts a route from the microscopic Hub- 
bard model on the square lattice to an effective lower 
energy Plaquette Boson Fermion Model (PBFM) at low 
hole doping. We apply the Contractor Renormalization 
(CORE) method of Morningstar and WeinsteinO to the 
plaquettized lattice (see Fig.^j) in a one step transforma- 
tion. 

We find that the bosonic part of the effective Hamilto- 
nian is-,glosely related to the projected SO(5) (pSO(5)) 
theorytffl; a theory of four bosons: a hole pair and a (an- 
tiferro)magnon triplet. The pSO(5) model describes the 
competition between antiferromagnetism and supercon- 
ductivity in a quantum mechanical framework. Its mean 
field theory yields some broad features of the cuprate 
phase diagram at low doping. At low temperature, 
hole pairs are governed by a phase fluctations actior 
which explains the (non BCS) proportionality between 
superfluid density, transition temperature and hole con- 
centration. In the superconducting phase, the magnons 
are massive and give rise to an antiferromagnetic res- 



onance in neutron scattering. These massive magnons 
were argued to-produce a resistance peak series in Joseph- 
son junctionals. 
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FIG. 1. Local bosons and fermions on the plaquette 
lattice. The singlet RVB vacua are depicted as solid squares. 
Holes are depicted by circles. The triplets, single hole and hole 
pairs Hubbard eigenstates define the degrees of freedom of 
the effective Plaquette Boson-Fermion Model. Interplaquette 
couplings are computed using Contractor Renormalization. 



Nevertheless without a microscopic foundation, the fun- 
damental "mechanism" problem remains: What creates 
and holds together d-wave hole pairs in the presence of 
local repulsive interactions, without the benefit of re- 
tardation and phonons? Even assuming that hole pairs 
move coherently, what is their hopping rate, and is it 
of the same order as the Heisenberg exchange energy as 
assumed by the pSO(5) theory? 

Here we address these questions and afford the 
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pSO(5)theory and its phase diagram a microscopic foun- 
dation. We also include the hole fermions which provide 
gapless (nodal) excitations in the square lattice. Their 
bandstructure is obtained from previously published nu- 
merical results, and their coupling to the bosons is esti- 
mated by symmetry and microscopic considerations. 
The resulting Plaquette Boson Fermion Model (PBFM) 
describes two coupled charged systems (i) Bole pair 
bosons, which Bose condense below T c and induce a prox- 
imity gap on the holes, (ii) The hole fermions, which 
occupy small Fermi pockets around (±7r/2, ±7t/2) and 
have a large van-Bove peak in density of states near 
(±7r, 0), (0, ±7r), the "antinodal" points. 
We discuss the thermodynamics of the coupled system, 
with the constraint on the total doping density. Previ- 
ously proposed boson-fermion models!!! differ from the 
PBFM by their Bilbert space (e.g. by having a large 
Fermi surface, and counting occupations from the elec- 
tron vacuum). 

For the PBFM, in the weak coupling approximation, 
some straightforward experimental implications are ob- 
tained: 

1. Hole spectral weight in Luttinger theorem- violating 
momenta (outside the "large" electron Fermi sur- 
face), e.g. on the line (tt, 0) — * (tt, 7r)E3. This weight 
survives above T c and is associated with excited 
holes moving in the correlated RVB vacuum. 

2. Asymmetry in tunneling conductance. At low dop- 
ing, x « 1, particle- hole symmetry is expected 
to be violated i.e. the positive bias conductance 
(injection of electrons) is suppressed by a factor 
proportional to x, relative to the negative bias con- 
ductance (injection of holes). Such a trend indeed 
appears in tunneling datal!3. 

3. The pseudogap doping dependence. The pseudpgap 
energy in photoemmissionEjE-3 and tunnelinglj is 
at the van-Bove peak of antinodal fermions. The 
decrease of pseudogap with doping follows the in- 
crease in fermion chemical potential. Its derivative 
with respect to doping measures the combined hole 
fermions' and hole pair bosons' compressibilities. 

4. Nodal transverse velocity. The quasiparticles 
proximity gap near the nodal directions deter- 
mines their transverse velocity. This velocity 
can be measured by photoemmission and optical 
conductivity!!!! We expect it to vanish at T c , and 
to be proportional to the Bose condensat e orde r 
parameter. Thus it should scale as vj_ oc yjT c (x). 

5. Hole dependent entropy. At temperatures above 
the superconducting transition, hole pairs evapo- 
rate into hole fermions, because of the difference 
in theirj-density of states. The doping dependent 
entropyliU is dominated by the fermion contribu- 
tion. 



The paper is organized as follows: In Section II we intro- 
duce the eigenstates of the Bubbard model on a plaque- 
tte. The local bosons and fermions are defined as the cre- 
ation operators of these eigenstates. The physics learned 
from the four site problem is instructive: Undoped, the 
ground state is a local resonating valence bonds (or pro- 
jected (i-wave BCS) state. The four bosons create the 
lowest triplet and the hole pair singlet. There are two 
degenerate spin half plaquette fermion states with sym- 
metry (tt, 0) and (0, tt). 

It has long been appreciated that in the Bubbard model, 
two holes cannot bind on a dimer bond, but they can 
bind on a plaquette (and on larger clusters)Ea. The hole 
pair wave function has d x i_ y i symmetry for 7r/2 rota- 
tions. The next step is to compute their intcrplaquette 
hopping rate in order to see whether they can preserve 
their integrity on the infinite square lattice. A short dis- 
cussion is included about plaquette "vacuum crossing" , 
which occurs at large chemical potentials and may be as- 
sociated with a transition from underdoped to the over- 
doped regime. 

In Section III the Contractor Renormalization (CORE) 
method is reviewed. The method requires exact diago- 
nalization of multi-plaquette clusters, in principle up to 
infinite range. Of course, the method is useful only if it 
converges rapidly in a feasible range of interactions. We 
have tested the convergence of the low spectrum for the 
Bubbard models on open ladders, with satisfying results. 
These tests confirm our belief that the convergence de- 
pends on a short boson coherence length, of order one 
plaquette size. This is very encouraging for the useful 
application of CORE to our problem, since the experi- 
mental superconducting coherence length of cuprates also 
appears to be particularly short in the underdoped sys- 
tems. We discuss the artifacts of the formal translational 
symmetry breaking within CORE. In Appendix^, we use 
the tight binding model as a pedagogical example of how 
longer range interactions of CORE serve to restore an 
unphysically broken symmetry. 

In Section IV the Plaquette Boson-Fermion model is de- 
rived. We discuss the hole pairs integrity, as evidenced 
from the numerical results, and how it is related to the 
sizeable pair hopping energy. The pair kinetic energy 
is crucial in stabilizing superconductivity. The full four 
boson hamiltonian is given in Appendix |aJ. The hole 
fermions band structure and interactions with the bosons 
are added. The thermodynamics of the weakly coupled 
PBFM yields a relation between the pseudogap energy, 
the bosons and fermions compressibilities, and the evap- 
oration of hole pairs into fermions at higher temperature. 
We conclude with a summary and a discussion of future 
directions in Section V. 



II. PLAQUETTE STATES 

We study the Bubbard model 
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where c- electron creation and number operators 

at site i on the square lattice. We will occasionally refer 
to its Gutzwiller projected version, the t-J model: 

n" = -tv ]T (is- s + "•'•) v • J X>-^ • J'- 

(ij),s {ij) 

(2) 



V is the projector of doubly occupied states, and J 
U 2 /U at large U/t. J' is the term of ordex, J that i 



eludes next nearest neighbor hole hoppingEB. For one 
electron per site (half filling), the short range antiferro- 
magnetic correlations are apparent when diagonalizing 
(|l|) and (||) on two sites. The dimer states, were used 
to construct effective models on the ladderEJ'Eil and for 
spin-Peierls phases on the square latticeEJ. The projected 
SO (5) model was defined on a. ladder using empty dimer 
states as the hole pair bosonsE However, there is no hole 
pair binding for the Hubbard model on a dimer. Naively, 
this suggests that pairs could readily disintegrate into sin- 
gle holes once inter-dimer hopping is turned on. More- 
over, if one wishes to capture (i-wave symmetry in the 
hole pair wavefunction, the basic unit block must pos- 
sess at least four-fold rotational symmetry. 
The smallest such block that can cover the square lattice 
is the four site plaquette. It is a trivial task to diago- 
nalize the Hubbard model on a plaquette and obtain its 
spectrum and wavefunctions. 
The spectrum is depicted in Fig. 0. 
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FIG. 2. Lowest spectrum of the Hubbard model on 
a plaquette. Eigenstates are labelled by total spin S and 
plaquette momentum Qx,Q y = 0, tt. The shaded area is over 
all high energy truncated states. The vacuum is denned as 
| fi ) , and the second quantized operators connect the vacuum 
to the lowest eigenstates as shown. 



lows (i) Real space (RS) description using holes, dimer 
singlets and dimer triplets as depicted in Fig. || (ii) 
Plaquette momenta (PM) representations using Q = 
(Qx,Qy), Qa = 0,7r, the four points on the plaquette 
Brillouin zone. The plaquette electron operator of pla- 
quette i is given by 



c f =- 



E 

rf—0,x,y,x-\-y 



p iQ-vJ 



(3) 



It is instructive to examine the plaquette eigenstates and 
energies in some detail before proceeding to couple them. 



A. The vacuum 

The ground state of the 4-site Hubbard model at half 
filling (n e — 4) is called | fi ) . In the PM representation 
it can be described by (suppressing the plaquette index), 



jfi) 



=(4 



o)TV,o)l 



C (0,7r)T C (0,7r)l) C (0,0)T C (0,0)l 



|o>. 

(4) 



Z is the wavefunction normalization factor. | fi ) is a d- 
wave BCS state, where doubly occupied states are sup- 
pressed by a partial Gutzwiller projection V{U jt). (At 
large U, V becomes a full projection). 
In the RS representation, see FigJ|, | fi ) is depicted as 
the resonating valence bonds (RVB) ground state of the 
Heisenberg model plus small contributions from doubly 
occupied sites. In the two dimer basis, | fi ) contains a 
large contribution from a triplet pair. 
The product state | fi ) = Yl^ laq \ fi ) i , is our vacuum 
state for the full lattice, upon which Fock states can be 
constructed using second quantized boson and fcrmion 
creation operators. 
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FIG. 3. Real space representation of plaquette 
bosons. Dominant spin and charge correlations in the pla- 
quette bosons wavefunctions. Bold lines represent singlet 
dimers t;J._, — Ulj, and double lines represent the triplet 
T4j + litji TiTi ± Ulj, where i and j are on sublattices A 
and B respectively. Holes are depicted by open circles. 



Since it is cumbersome to write the full wavefunctions ex- 
plicitly, we represent their dominant correlations as fol- 
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B. Magnon triplet 

The magnons arc defined by the lowest triplet of S = 1 
states. In PM representation they are 



4 |n) =-^Y, c Qs<>c Q+(7r ^ s ,\n) 



Qs 



a = x,y,z, 



(5) 



where o~ a are Pauli matrices. These (antiferro-) magnons 
have plaquette momentum Q = (tt, it). Their excitation 
energy is close to the superexchange energy J w 4t 2 /U. 
An antiferromagnetic state can be constructed by a prod- 
uct of plaquette coherent states 



plaq 



(6) 



where |m| = 1. This state supports a finite staggered 
moment 



I 

N 



= ■\/3/STn a cos 9 sin 6 

< 0.306. (7) 



Note that the maximal magnetization per site supported 
by \I/ a / m is less than the classical value of 0.5, since it 
does not contain higher spin states up to S = 2. 



C. Single hole fermions 

The ground states for a single hole (n e = 3) are two 
degenerate doublets described by plaquette momenta 
Q=(0,7r),(7r,0): 



v 



mi 



=T-L. 



(8) 



state has vanishing hole spectral weight, that is to say 
for all values of U/t 



Z (ir,ir) — \{Q\f{-K,-K)s c {ir,-K)s I n) | 2 — 0. 



(11) 



Since these states couple by interplaquette hopping to 
the lower doublet, this produces an asymmetry of the 
quasiparticle weight between momenta close to (0, 0) and 
(tt, tt) . This asymmetry may explain the difficulty in ob- 
serving "shadow bands" , i.e. quasiparticles on the Fermi 
pockets surfaces closer to (tt, np3. 

It is interesting to note that the two-fold degeneracy of 
the fermion doublets is a property of the plaquette. The 
four site Hubbard and t-J Hamiltonians happen-io com- 
mute with the plaquette d-density wave operatorE3 



D 



0) s C(0,tt)s ~ c (0, w ) s C (t,0)s)'^- 



(12) 



D connects between the doublet pairs (tt,0) <-> (0,7r), 
and (0, 0) <-> (71", tt). Thus a possible ground state of one 
hole is the current carrying state 
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(0,7T) 



mi 



(13) 



which is a staggered flux (or d-density wave) state. For 
a single hole in 4 x 4 periodic lattices, t his sta te does not 
seem to be the lowest energy (see Section IV D). However, 
a large susceptibility for such currents is expected since 
the hole dispersion has a valley between the antinodal 
points (tt, 0) and (0,7r), which is weakly dispersive and 
contains a large admixture of the two plaquette fermion 
states. It is thus conceivable that the staggered flux com- 
bination of the plaquette fermions would be selected in a 
vortex core or near the sample edge. 



D. Hole pair boson 



where . . . represent higher order electron operators. The 
hole fermion Bloch state can be constructed as 



t 

k+Qs 



plaq 
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e »k. Xl/ t 



Qis 



m: 



(9) 



For a lattice of disconnected plaquettes, /J Q y creates an 
eigenstate with a photoemmission spectral weight given 
by 



\(®>\f(n,0) S C(Tr,0)s I n) | 2 = Zfrfi), 



(10) 



where e.g. for the t-J model, 1/4 < Z/ n ,o) < 1/2 is a 
function of tj J. This weight is further renormalized by 
interplaquette couplings in the effective Hamiltonian. 
Incidentally, there is another degenerate pair of dou- 
blets at higher energy (of order J) at momenta Q = 
(0,0), (tt, tt). It turns out that by symmetry, the (tt, tt) 



The ground state of two holes (n e — 2) is described by 



--Vc 



(0,0)^(0,0)1 



|0) 



^<%c iT c jU + ... |fi), (14) 



where d+j is +1 (-1) on vertical (horizontal) bonds, and 
. . . are higher order [7/t-dependent operators. Thus, b> 
creates a pair with internal d-wave symmetry with re- 
spect to the vacuum. For the relevant range of U/t, the 
state normalization is 1/3 < Z' b < 2/3. The important 
energy to note is the pair binding energy defined as 



A b =E(0)+E(2)-2E(l) 



(15) 



where E(Nh) is the ground state of Nh holes. Ah is 
depicted in Figure 0. In the range U/t 6 (0,5), it is 
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bounded by — 0.04< < Af, < 0. It has been well ap- 
preciated that the Hubbard, t-J and even CuC>2 mod- 
els have pair binding in finite clusters starting with one 
plaquetteEa. In larger clusters, such as the 4x4 lattice, 
pair binding is seen for up to 6 holesa (three hole pairs) 
for U/t < 20. This does not yet explain the integrity of 
pair correlations on the infinite lattice since the electron 
hopping energy t is much larger than the pair binding en- 
ergy. In Section IV B and Appendix |X], we show numer- 
ical evidence that plaquette pairs survive disintegration 
into fcrmions. 
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FIG. 4. Pair binding energy on a plaquette. A& of Eq. 
( [til ) calculated for the Hubbard model for different interaction 
strengths. For Aj < the hole pair is more stable than two 
single holes on a disconnected plaquette lattice. 



A d-wave superconducting state can be written as the 
coherent state 



plaq 



<Jf u 



-scF 



U(cos0 + sin0e i¥, &t) |fi) , (16) 

i 

with the superconductor order parameter 



order parameter. It should have larger matrix elements 
than the customary d-wave pairing operator dijC^Cji. 



E. Underdoped to Overdoped Transition 

Throughout this paper we restrict ourselves to low dop- 
ing, i.e. a small number of hole pairs and hole fermions 
per plaquette. Nevertheless, the plaquette states basis 
leads us to expect an interesting transition at higher hole 
doping for the following reason. 

When the chemical potential is large enough to bring the 
hole pair (2 electron) state to be lower than the 4 elec- 
tron vacuum, a vacuum crossing takes place. On a single 
plaquette, the vacuum crossing is when the levels of zero 
and one boson intersect. Since hole pairs are somewhat 
larger in size than a single plaquette, the lattice vacuum 
crossing should take place at somewhat less than x = 0.25 
holes per square lattice site. 

Once the two hole state turns into the new vacuum | f2 ) ', 
all excitations are defined with respect to it using differ- 
ent boson and fermion creation operators. For example, 
the old RVB vacuum becomes a Cooper pair excitation 
above the new vacuum. 



in: 



(Vy due 



(18) 



It is plausible that the vacuum crossing is a true quantum 
phase transition, and not merely a mathematical artifact 
of using different plaquette bases to construct the same 
ground state. A candidate for such a phase transition is 
the restoration of square lattice symmetry, if this symme- 
try is truly broken by plaquetti zation in the underdoped 
regime as mentioned in Section III B . 
The overdoped side is far from half filling, where effects 
of the Gutzwiller projection are small. That is to say, 
the eigenstates can be approximated by applying elec- 
tron operators to the electron vacuum. With interpla- 
quette hybridization of the two electron plaquette vacua, 
the ground state can be adiabatically connected to the 
quarter filled electron Fermi surface. In the absence of su- 
perconductivity it will exhibit a large (Luttinger-theorem 
obeying) Fermi surface. 



{•$>\dijC, 



Z' b e %{p sin^ cos 9. 



(17) 



It is worthwhile to reemphasize the following point which 
has important implications in interpreting both exper- 
iments and numerics of Hubbard like models. The 
fermions and hole pairs have charges +e and +2e respec- 
tively. Their number is counted from the correlated half 
filled (RVB) vacuum. The operators / and b should not 
be confused with the electron and Cooper pair operators 
c\ and ejej respectively, whose numbers are counted up 
from the electron vacuum. In numerical calculations, it 
is preferable to use the operator b', as defined by the 
Hubbard plaquette eigenstates, as the superconducting 



III. CONSTRUCTION OF THE EFFECTIVE 
HAMILTONIAN 

Having described the low lying plaquette states, we are 
faced with the challenge of constructing an effective 
Hamiltonian for the full lattice. Motivated by the pair 
binding on a plaquette, one might initially wish to com- 
pute the effective hopping of a hole pair between pla- 
quettes using second order perturbation theory in the 
interplaquette hopping t' . 

This naive approach yields pair hopping of order J c oc 
The pcrturbativc expansion is controlled by t'/A&. 
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This suggests failure of perturbation theory for t' = t > > 
Afc. Indeed, by looking at the exact spectrum of two 
connected plaquettes, we find that second order pertur- 
bation fails in a sizable domain of if < t. Does this imply 
dissociation of the local bosonic correlations in the trans- 
lationally invariant lattice? 
We argue no. 

There is convincing evidence from various numerical ap- 
proaches, that two holes remain on the same plaquette 
in -\/26 x V26 latticesn and on 8 x 6 t-J laddersa. How- 
ever, in order to understand how many pairs behave on 
the infinite lattice, we must determine the pair hopping 
energy, and derive their effective Hamiltonian. 
A suitable approach for this task is provided by the Con- 
tractor Renormalization (CORE) methodO described be- 
low. The small parameter of CORE is the ratio of the 
hole pair separation, i.e. coherence length to the range 
of the effective interactions. 



Representing Ji ren in the real space plaquette basis 
| Oil, ... , ckjv ) , defines the (reducible) interplaquette cou- 
plings and interactions. 



H, 



ren 



'(1,2,3) 




M states 



^=\> 



FIG. 5. The reduced Hilbert space of a plaquette 
cluster used in CORE. The calculation of Hl^s requires 
diagonalization of the Hubbard model on the cluster high- 
lighted by the shaded region. }{ Ten reproduces the exact 
spectrum within the reduced Hilbert space. 



A. Contractor Renormalization Algorithm 

Given a microscopic Hamiltonian 7i on the square lat- 
tice we choose a plaquette covering and proceed by the 
following steps: 

Step 1: Defining the reduced Hilbert space. We diagonal- 
ize H on a single plaquette and truncate all states above 
a chosen cutoff energy. This leaves us with the lowest 
M states { | a) }f f . The reduced lattice Hilbert space is 
spanned by tensor products of retained plaquette states 
| oil, ... , ctN ) ■ A case in point is the Hubbard model 
spectrum, which for the half filled case has 70 states. We 
truncate 66 states and keep the ground state and lowest 
triplet, i.e. M = 4. Thus, the Hilbert space is consider- 
ably reduced at the first step. 

Step 2: The Renormalized Hamiltonian of a cluster. The 
reduced Hilbert space on a given connected cluster of N 
plaquettes is of dimension M. — M N . See Fig. |5| for an 
illustration. We diagonalize TL on the cluster and ob- 
tain the lowest M. eigenstates and energies: (\n) ,e„), 
n = The wavefunctions |n) are projected 

on the reduced Hilbert space and their components in 
the plaquette basis | a%, . . . ) are obtained. The pro- 
jected states tp n are then Gramm-Schmidt orthonormal- 
ized, starting from the ground state upward. 

\$») = Y ( IV>n) - £ \M ^rnHn) ) , (19) 

where Z n is the normalization. The renormalized Hamil- 
tonian is defined as 



M 



(20) 



which ensures that it reproduces the lowest M. eigenen- 
ergies exactly. 



Step 3: Cluster expansion. We define connected N point 
interactions as: 



jrren 

U {i 1 ,...,i N ) 



(ii,—,i' N ) 



(21) 



where the sum is over connected subclusters of 
. . . , In). The full lattice effective Hamiltonian can 
be expanded as the sum 



ff 



(ijk) 



(22) 



hi is simply a reduced single plaquette hamiltonian. hij 
contains nearest neighbor couplings and corrections to 
the on-site terms hi. hijk contains three site couplings 
and so on. hi t ^_ t i N will henceforth be called range-N 
interaction. We expect on physical grounds that for a 
proper choice of a truncated basis, range-N interactions 
will decay rapidly with N. This expectation needs to be 
verified on a case by case basis. 

Morningstac, and Weinstein, by retaining up to range-3 
interactionsa, demonstrated that the CORE renormal- 
ization group flow, obtains an excellent value for the 
ground state energy of the spin-| Heisenberg chain. This 
is encouraging, since the spin half chain has long range, 
power-laai decaying spin correlations. Pieckarewicz and 
ShepardEj tested CORE for the 12 site spin-i Heisenberg 
ladder. They got better than 1% accuracy for all lowest 
64 states using a plaquette basis keeping only up to range 
two interactions. 

In general, there is no apriori quantitative estimation of 
the truncation error. Nevertheless, if it decays rapidly 
with interaction range, we deduce that there is a short 
coherence length related to our local degrees of freedom, 
e.g. in our case the hole pair bosons and the triplets 
(bound states of two spinons). 
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1. CORE of Wave Function Correlations 

The CORE process is designed to reproduce the low ly- 
ing spectrum, while the wave functions may be signif- 
icantly distorted during the truncation of the Hilbert 
space. Does this hinder calculation of correlation func- 
tions within this scheme? The answer depends on the 
operators whose correlations we wish to calculate. 
The correlations are calculated by adding external source 
terms to the original Hamiltonian, 



H[r)] = H[0}+J2vidi 



(23) 



where Oi is a microscopic linear perturbation operator 
whose correlations we wish to determine. 
We apply the CORE cluster expansion of Eq. (|2^ ) to 
the perturbed Hamiltonian and obtain H e ff [77] which is 
expanded to linear order in rji, 

KeffM = E { h i + vA?) + E (>• + nid$ 

3 (jk) 
(jkl) 

^K ff + J2vA ff + o(v 2 ) 



A rapid decay of hi 1: ... ln beyond a short truncation range 
is essential for the feasibility of the CORE scheme for 
the spectrum. Similarly, to calculate the correlations of 

0, we require a rapid decay of with n, which would 
allow us to calculate 0\ by small clusters diagonaliza- 
tions. We have previously argued that a small trunca- 
tion error results from a short coherence length £, for the 
coarse grained degrees of freedom. For example, the size 
of the hole pairs in the slightly doped Hubbard model. 
Thus, the operator O should be chosen to have large ma- 
trix elements within the reduced Hilbert space, and the 
multi-site operators t should decay rapidly beyond 
the range of £. 

In other words, if the truncated wave functions retain the 
relevant local operator content, the cluster expansion for 
the operators converges rapidly, and the effective Hamil- 
tonian can reproduce the long wavelngth correlations cor- 
rectly. In this paper, the truncated plaquette states, for 
example, contain d-wave hole pairs on plaquettes. If 
these hole pairs turn out to be tightly bound in the ex- 
act eigenstates of the full lattice, their creation operator 
has small multi-site corrections in the renormalizcd basis, 

1. e. it has a rapidly decaying cluster expansion. In this 
case, long range d-wave pair correlations are well repre- 
sented (up to an onsite renormalization factor), by the 
boson-boson correlations of the Four Bosons model. 



°r // = o| 1) +E^ ) +E°S + 



(24) 



jk 



0\ represents the linear perturbation Oi in the trun- 
cated Hilbert space, 

M 



6f f =^2(n\di\m)\^ n ) (i> m \ 



(25) 



The two-point dynamical correlation function at low tem- 
perature T « e m ax is given by the Lehmann represen- 
tation 

Sij(uj) = 

2yj- e 

-— }^e"^(n\Oi\m)(m\Oj\n)5(u> + e„ - e m ) 



|^ E e-^^ n \dt f, \^ m )(^ m \df f \^) 



trunc 



ZN 



x 5(w + e„ - e m ) 



(26) 



where Z is the partition function, and the second sum 
is evaluated in the truncated Hilbert space, using the 
cluster expansion (j24|) for the matrix elements. 
Thus for the full lattice cluster, O e ^ recovers the ex- 
act correlations of the true low energy eigenstates. By 
( p4| ) the linearized cluster expansion involves multi-site 
operators 



O 



(„.) 



dhi lt ... ln /dr)i 



(27) 



B. Lattice Translational Symmetry 

The CORE algorithm formally requires explicit break- 
ing of lattice translational symmetry at the first step. 
The plaquette lattice vacuum breaks lattice translational 
symmetry as follows: each plaquette vacuum contains a 
triplet pair contribution, but the product state does not 
contain interplaquette triplets, and hence differs from the 
state translated by one square lattice spacing. In order to 
restore the lattice symmetry, interplaquette triplet pair 
correlations can be reintroduced by triplet pair creation 
operators in the effective Hamiltonian. 
As an illustration, let us consider the plaquette vacuum 
of the two leg ladder, in Fig.^| which can be written in 
the form 



1 



|0) = Vexpl-^ 



t 



2ia u 2i+l,a 



ni°>*> 



(28) 



where |0)j is a singlet on rung i, and t| a |0)i is a rung 
triplet. V projects out multiple occupation of triplets. 
The translational invariant RVB state in FigJ^, can be 
constructed from | O ) by applying the operator 



|RVB}=7>exp ^E4-iAJ 



(29) 



In the triplet bosons representation of the Heiscnberg ex- 
change there are anomalous inter-plaquette terms t\t\ + 
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FIG. 6. Restoration of translational symmetry on 
the ladder. fi ) is the plaquette lattice vacuum which 
breaks two fold lattice translational symmetry. |RVB) is the 
dimer Resonating Valence Bonds state which has translational 
symmetry, and is related to | Q ) by an exponential of triplet 
pairs, see Eq. (ps|) ■ 



titj. In the mean field theory the Bogoliubov transforma- 
tion introduces a partial "symmetry restoring" exponen- 
tial operator, which resembles (|29|). However symmetry 
cannot be fully restored to the wavefunctions in the re- 
duced Hilbert space because of the elimination of higher 
spin states. 
Is this a problem? 

Well, it depends on what one is interested in. CORE 
is constructed to obtain accurate effective interactions. 
Unphysical symmetry breaking effects can be introduced 
by truncating the longer range interactions, as explained 
by a toy model in Appendix |b| Therefore it is hard to 
rule out a physical "plaquettization" of the true ground 
state. Incidentally, such a fourfold discrete symmetry 
breaking is consistent with Berry phase argumentsEj for 
the spin liquid phase of spin half Heisenberg models. 
The symmetry breaking, appears as minigaps near the 
edges of the plaquette lattice Brillouin zone (PLBZ) 
n x ,Ky € (— 7r/2, 7r/2). In Appendix |b], we see how the 
minigaps of the tight binding model, decrease as longer 
range interactions are included. For the triplet and hole 
pair bosons, minigaps do not matter much since their low 
energy states are around (0, 0) and (ir, ir) respectively; 
the farthest possible from the PLBZ edges. 
On the other hand, low energy fermions happen to be 
centered around the PLBZ corner (±7r/2, ±7t/2), where 
the effects of plaquette symmetry breaking on the spec- 
trum are large. Although by rotational symmetry, the 
two bands which contain the (tt, 0) and (0, tt) states are 
degenerate at the PBLZ corner, the other two bands have 
minigaps. These would distort the elliptical shape of the 
Fermi pockets, an effect which if it exists, could be de- 
tected by angular resolved photoemmission. 



IV. THE PLAQUETTE BOSON-FERMION 
MODEL 

We first start with the bosons, and compute their inter- 
plaquette couplings and interactions using CORE. Later 
we introduce the hole fermions, whose parameters are 
taken from published numerical data on large clusters, 
and estimate their coupling to the bosons using symme- 
try arguments. Finally we discuss the properties of the 
combined Hamiltonian. 



A. Computing Boson Interactions 

For the purpose of this paper, we have limited the CORE 
calculations to range-2 boson interactions, while project- 
ing out the fermion states. This required a modest nu- 
merical diagonalization effort of the Hubbard model on 
up to 8 site clusters. The resulting range-2 Four Boson 
model can be separated into bilinear and quartic (inter- 
action) terms: 

n ib = n b [b] + u l [t] + n mt [b, t] (30) 

where the bosons obey local hard core constraints 

&I*>i + X>L*ai<l (31) 
a 

The bilinear energy terms are 

n b = (e b - 2/*) b \ b r ~ J bJ2 { b h + H - c -) 

= e t -yE (4*<« + H - c -) 

~ J f 5>L^+H.c.). (32) 

a{ij) 

In Fig. [7] we compare the magnitudes of the magnon hop- 
pings Jt, Ju and the hole pair hopping for a range of 
U/t. First, we observe that J t « Ju ~ 0.6J, i.e. the 
magnon terms have similar form as those previously ob- 
tained for the HeisenbeEft model in the bond operatorEJ, 
and plaquette operatorEO representations. Second, the 
region of intersection near U/t — 8, is close to the pro- 
jected SO (5) symmetry point. We emphasize that al- 
though there is no quantum SO (5) symmetry in H 4b , 
there is an approximate equality of the bosons hopping 
energy scales. r-This equality which was assumed in the 
pSO(5) theorya, previously appealed to phenomenologi- 
cal considerations. Here, the equality emerges in a phys- 
ically interesting regime of the Hubbard model and has 
important consequences on the phase diagram as shown 
below. 

Ti mt includes nearest neighbor triplet-triplet, pair-pair, 
and pair-triplet interactions. In Appendix TC 4b with 
all its terms is displayed in its full glory, and a table of its 
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FIG. 7. Boson hopping energies versus Hubbard U. 

Jt, and J t t are the magnon's normal and anomalous hopping 
energies. Jt is the hole pair hopping energy. The intersection 
region near U — 8 is close to the projected SO(5) symmetry 
point. All energies are in units of t. 



computed coupling constants is provided for the values of 
U/t = 3, 6, 8, 10. We also compute the truncation error of 
discarding range-3 terms. This is done by comparing the 
Hubbard model at U/t = 6, with and 2 holes on 12 sites, 
to that of corresponding range- 2 H 4b . Relative shifts 
of less than 1% in the ground state and first excitation 
energies, correspond to a very small truncation error. 



B. Mechanism of Superconductivity 

There are two important effects which together can lead 
to supercondutivity: (i) Pairing, and (ii) Bose conden- 
stion of the pairs. An important energy scale for both 
effects is the pair hopping rate Jf,. 

The small range-3 truncation error was found at a large 
interaction U/t — 6, where there is actually no pair bind- 
ing on a single plaquette (see Fig. ||). The convergence 
of effective interactions, implies short boson coherence 
lengths £,t,£,b- £t is the distance between spinons (local- 
ized spin half configurations) which comprise a magnon. 
£b is the hole pairing distance. Both coherence lengths 
appear to be of the order of one plaquette size. This con- 
clusion is supported by numerical observation of short 
distance (latticejfpjnstant) correlations between two holes 
on large latticeso-B. It is interesting that the short pair 
coherence length is dynamically generated in the Hub- 
bard model, even for U/t > 4.5 where the pair binding 
energy on an isolated plaquette is positive. 
Why is £ so short? 



There are two effects which bind pairs: a classical mag- 
netic energy from minimizing the number of broken 
Heisenberg bonds, and a quantum kinematic pairing for 
holes moving on two sublattices of a quantum disordered 
antifcrromagnet. The first effect is supported by finding 
pair binding on a single plaquette. However, this energy 
also favors clumping many holes together. The quan- 
tum pairing effect was proposed by Weigmann, Lee and 
Wen, who integrated out spin fluctuations in a quantum 
disordered phase, to induce a long range electrodyr, 
cal attraction between holes on opposite sublatticescal 
The kinematic effect produces pairing rather than phase 
separation, and is robust against additional short range 
repulsion. It also can explain-pair binding on large clus- 
ters in a regime of U/t > 4.5L3. 
Bose condensation. 

The relative large hopping in the pair kinetic energy 
— Jb^2u j\ b\bj is crucial for understanding the cuprate 
phase diagram. 

(i) The pair kinetic energy competes effectively with 
the antiferromagnetic order. While uncorrelated single 
fermion kinetic energy is not inhibited by the presence of 
long range antiferromagnetic order (in fact it strengthens 
it by a Nagaoka-like mechanism), the pair kinetic energy 
is substantially lower in a background of short range sin- 
glet correlations. This effect was clearly demonstrated 
in variational Monte Carlo studies of pair kinetic energy 
in doped RVB wavefunctionsEl, and is also a property of 
the variational treatment of the four boson model. 

The destruction of antiferromagnetic order into a quan- 
tum spin liquid with massive triplets, also helps in the 
kinematical pairing process as discussed above. 

(ii) Having destroyed antiferromagnetic order, the pair 
kinetic energy competes with charge localization due to 
disorder, or solidification (charge density wave), and with 
disintegration into unbound hole fermions. 

(iii) A large Jb stabilizes a superconducting phase at finite 
temperatures. It determines the superfluid density p s ^, 
2J C |(6)| 2 , and the phase ordering transition temperatures 
T c ps p s . 



C. Four Boson Mean Field Theory 

The mean field theory is separated into two parts: (i) 
Calculation of the order parameters as a function of dop- 
ing, using variational coherent states, (ii) Determina- 
tion of magnon resonance energy from a soft interaction 
version. The results are qualitatively similar to the pro- 
jected SO (5) phase diagramtffill. 

Here, we choose U — 8t, and evaluate the energy of the 
full boson Hamiltonian (A3,A4) in the variational coher- 
ent states ip afm (8) and ^ 



^6) of Eqs. (|6|) and (|16|) re- 
spectively. These states represent the antiferromagnetic 
and superconducting phase. The critical chemical poten- 
tial He, where the ground state energies cross, and quan- 
tum fluctuation angle 9(fi) arc determined by minimizing 
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the energy. The spin stiffness and superfluid density are 
given respectively by 



Paf = 2J t (t) 2 , 
Psc = 2J b (b} 2 , 



(33) 



where we use Eqs.(0,0) for the magnon and hole pair 
expectation values. These coefficients, which determine 
the transition temperatures, as well as the doping concen- 
tration x are plotted as a function of chemical potential 
in Fig.|[ We emphasize that the results should not be 
quantitatively compared to experiment, since they are 
variational approximations to a simple model, and ne- 
glect effects of low energy hole fermions. 
The variational theory yields a first order transition be- 
tween zero doping and x c = x(p c ) 0.125, where the 
staggered magnetization abruptly vanishes and the su- 
perfluid density jumps to a finite value. 
For charged holes, this first order transition (phase sepa- 
ration), is forbidden by long range Coulomb interactions. 
Instead one expects high compressibility, incommensu- 
rate mixed phases and stripesE 2 ] in the intermediate dop- 
ing regime x € (0,x c ). 

Even a weak disorder potential is very efficient in b, 
ing the intermediate phase into "quantum melts" 
i.e. puddles of superconductor inside antiferromagnetic 
domains. 

Above x c , the superfluid density increases with dop- 
ing, in agreement with London penetration depth mea- 
surerhentsu, The overdoped regime is beyond the ex- 
pected vacuum crossing point (see Section HE). 
The magnon dispersion in the superconductiruz-chase is 
obtained by decoupling a soft core interactionOcil, 



n 



int 



(34) 



where W is fitted to yield the order parameter magni- 
tudes calculated variationally. 

In the superconductor, the magnons acquire a gap at the 
antiferromagnetic resonance u res which increases with 
doping as 



uj res = 2\J (p - p c ){p - p c + 2J t ) oc ^x - x c . (35) 

This dependence, as plotted in FigJ|, is qualitatively con- 
sistent with inelastic neutron scattering dataE3. 



D. Fermion Hamiltonian 

In the previous section we have computed the bosonic in- 
teractions of ( |30| ) from the Hubbard model using CORE. 
In that computation, we have eliminated the fermion 
(single hole) states, we expect however that for the two 
dimensional square lattice, low energy fermion excita- 
tions are important. While the fermion holes short range 
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FIG. 8. Variational Solution of the Four Boson 
model. Results correspond to Hubbard interaction strength 
U/t = 8. the antiferomagnetic resonance energy, paf 

is the spin stiffness in the antiferromagnetic phase, and psc 
is the superfluid density in the superconducting phase, p — p c 
is the chemical potential difference from the first order tran- 
sition at p c , and x is the hole density. The estimated vacuum 
crossing point is discussed in Section |ll E|. 



effects on the boson couplings were included in the range- 
2 CORE calculations, their long wavelength excitations, 
require diagonalizing larger clusters which are beyond 
this paper's computational scope. We therefore resort 
to including the hole fermions dispersion "by hand" i.e. 
use the single hole band structure computed previously 
for large clusters. We then estimate their interactions 
with the bosons. 

It is important to emphasize that the definition of the 
hole pair bosons and the hole fermions is simply a matter 
of separation: two hole fermions are on different plaque- 
ttes. When they hop into the same plaquette they turn 
into one boson via the Andreev coupling defined below. 
For the relevant range of U ft the numerically determined 
band structures for the single hole can be fit by two hop- 
ping energies 



ks 



t' (cos(k x a) + cos(fc y a)) 2 

+ t" (cos(k x a) — cos(fc y a)) 2 



(36) 



k runs over the square lattice Brillouin zone. See Fig. |. 
The values t' m J, and t" w 0.1J are taken from the 
numerical Quantum. Monte Carlo data for the t-J model 
on a 24 x 24 latticeE3, find that for the physically relevant 
range of J £ [0.4t, 0.6£], the dispersion values are t' ~ 
0.7 J, t" ~ O.lt'. 

The magnitude of t' « J (rather than the bare value t) 
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This has important experimental implications: 



ant in ode 




FIG. 9. Hole fermions bandstructure. Contours of the 
single hole spectrum Eq.(^ ) modelled by fitting to published 
numerical results. For a dilute number of free holes, Fermi 
pockets will be created around the points (±7r/2, ±7t/2). A 
flat valley near the magnetic zone edge dominates the low 
energy spectrum. 



and the position of the minima on the magnetic zone edge 
(-7T, 0) — (0, tt) were explained by theories oOioles in the 
short range antiferromagnetic environments. The semi- 
classical theoryEa finds that holes are highly dressed local 
spin polarons, which effectively hop on one sublattice. 
From the CORE's perspective, the flat valley between 
(tt, 0) and (0, tt) is related to the original degeneracy be- 
tween the two lowest plaqucttc fermions. Thus we expect 
the wavefunctions of the fermions on the lattice to con- 
tain a large component of these two states. Consequences 
of this on the quasiparticle weight and possi ble s taggered 
orbital currents were mentioned in Section [IC. 



The holes have hard core interactions among themselves, 
and with the bosons. At low doping however it is still 
meaningful to describe their states by excitations about 
small Fermi pockets around (±7r/2, ±7r/2). 
The fermion density of states of ( pq ) is plotted in Fig. 0. 
We see a large peak at low energies (of order 4i" << 4r) 
from the saddlepoints near the antinodal points. These 
dominate the hole spectral function, and tunneling den- 
sity of states at the "pseudogap" A pg energy above the 
chemical potential. Within this framework, A pg does not 
describe the pairing correlation per se. (It only feels the 
change in boson density through changes in the common 
chemical potential). Even in the superconducting phase 
where hole pairs Bose condense, near antinodal points the 
Bogoliubov particle-hole admixture is small, and quasi- 
particles have a character of holes in the RVB vaccuum. 



1. Angular Resolved Photoemmission. The large 
Fermi surface of electrons, given by Hartree-Fock 
approximations, includes mostly the first magnetic 
Brillouin zone (the diamond connecting antinodal 
points). Luttinger's theorem for a Fermi liquid 
of electrons excludes any hole spectral weight out- 
side this area. In contrast, spectral weight of our 
fermions can be found anywhere outside the small 
Fermi pockets near (±7r/2, ±7r/2). Indeed, broad 
quasiparticle weight, above T c has been observed 
in photoemmission date, at momenta half way on 
the line (tt, 0) — > (tt, ir)£3. 

A direct evidence of small Fermi pockets would be 
sharp gapless quasiparticle modes on both sides of 
(tt/2,tt/2),. The "shadow" quasiparticles closer to 
(tt, tt) are harder to observe than the ones closer to 
(0,0), because of vanishing quasiparticle weight as 
discussed following Eq.(0). 



2. Tunneling conductance should exhibit an inherent 
asymmetry between injecting electrons (positive 
bias) and injecting holes (negative bias). The neg- 
ative bias peak at the pseudogap voltage is larger 
than the positive peak, since injecting electrons 
is suppressed by Hubbard interactions. In other 
words, electrons can only be injected into existing 
holes, whose density is of order x, at low doping 
the ratio of weights should scale with x. A re- 
view of (unsymmcti'iaed) tunneling data published 
by several groupsPtr 8 ! reveals such an asymmetry, 
although we have not seen yet a systematic study 
of its doping dependence in the literature. 



E. Boson-Fermion Couplings 

Couplings between bosons and fermions can be derived 
by microscopic considerations and symmetry. 



ELD 



j 



FIG. 10. Andreev coupling between hole pairs and 
fermions. The microscopic origin of the hole pair-hole 
fermion coupling is a simple unbinding process. Because of 
the hole pair d-wave symmetry, the coupling matrix element 
dij is odd under 7r/2 rotations on the lattice. 



Taking into account the c£-wave symmetry of the hole pair 
state yields an Andreev coupling (see Fig.|Io|): 



9b 

k,q 



t k+q/2 



k+ql 



H.c 



(37) 
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where — cos(k x ) — cos(k y ), and b q = ^? ag e iqXi &t is 
a Fourier component on the plaquette lattice. 
In the superconducting phase (b) ^ 0. This implies a 
proximity induced pairing of fermions in the small pock- 
ets, and an opening of a superconducting gap with the 
Bogoliubov dispersion 



± 



Ai c = 56d k (6). 



(38) 



To be consistent with the range-2 CORE method, we 
must not include close-by holes on nearest neighbor pla- 
quettes. These excitations were already taken into ac- 
count in the effective hole pairs hopping energy. The 
remainder Andreev coupling is therefore between second 
nearest neighbor hole fermions, with a coupling constant 
9b < 0.1 Jb, estimated from the magnitude of the range-3 
terms (see Appendix [A]). 

We emphasize that A sc is not the "usual" BCS gap, since 
it couples to hole fermions, not electrons. Through its de- 
pendence on the Bose condensate order parameter {b}rx, 
we can deduce the transverse quasiparticle velocity at the 
nodes v± = 9Ak/<9fcj_. v± — > at T c , and it should vary 
with doping as 



Vj_ oc 



'T c oc \fx. 



(39) 



At higher temperatures than T c , Ak vanishes and a 
broadened signature of the small Fermi surface emerges 
in the spectral function. In contrast, near antinodal 
points, Bogoliubov particle- hole mixing is negligble and 
spectral weight is due to hole fermions. This is consis- 
tent with photoemmission data which finds that above 
T c the gap, closes only in a small region around the nodal 
directionEl 

The Andreev coupling (|37]) couples the superconducting 
phase fluctuations to nodal quasiparticles. Similar inter- 
actions were used to calculate the temperature dependent 
London penetration lengthen. That calculation found the 
fermions to be more dominant at low temperatures than 
thermal phase fluctuations in destroying the superfluid 
density. The effects of this term on the fermions above 
T c , were recently arjpied to give rise to marginal Fermi 
liquid spectral peakscll 

Lastly, the fermion-magnon coupling is given by 

^ ^ — 9t + ^-m-q)/ks/k+q+7fs+m + H.C 

mskq 

(40) 

This singlet interaction term, flips fermion spins and 
scatters them with momentum (it, it) while emitting or 
absorbing magnons. It produces signatures of the , 
ferromagnetic resonance in the fermions self energj 
( pp| ) is similar to fermion-magnon terms which were con- 
sidered for predicting antiferromagnetiqajesonance signa- 
tures in tunneling and photoemmissionfca. 




FIG. 11. Hole fermions and hole pair bosons density 
of states. The fermion density of states is calculated for 
dispersions Eq. (^) for the normal state (dashed line), and 
superconducting state Eq.(^) (solid line). The low energy 
scale t" creates a large peak in the hole density of states. Hole 
pair bosons single particle density of states is approximated as 
a constant cor resp onding to the non interacting bilinear terms 
"'" ; " Eq. (]32j) . A sc and A pg refer to the superconducting 



of H" 



gap of eq. (38) and pseudogap of eq. (37) respectively. 



F. Boson-Fermion Thermodynamics 

The end result of the previous sections is a system 
of four bosons and a gas of hole fermions in ther- 
mochemical equilibrium, i.e. the charged bosons and 
fermions share a common chemical potential /i. Com- 
bining yields the complete Plaquette 
Boson-Fermion Hamiltonian: 



H PBFM = H 4b^ + n f^ +n bf +n 



If 



(41) 



In a uniform phase, the fermions and and hole pair bosons 
obey a global charge density constraint 



2n 6 (2^, T) + n/(/z, T) = x. 



(42) 



An important missing parameter, in the absence of a con- 
sistent calculation of the fermions bands, is the relative 
position of the lowest fermion and hole pair energies. 
Numerical evidence for 4x4 Hubbard clustersy show that 
for up to three hole pairs, there is a negative pair binding 
energy, i.e. the lowest fermion state at (tt/2,tt/2) is still 
above the boson condensate. However, at finite doping 
where superconductivity wins over antiferromagnetism, 
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the repulsively interacting bosons may have higher en- 
ergy than the bottom of the fermion bands. This will 
produce gapless nodal fermions in the superconductor. 
Here we shall assume that already at very low doping 
these energies match, and bosons and fermions coexist. 
The boson and fermion compressibilities are 



n b = drib/d(2fi), 
Kf = dnf/d^i. 



(43) 



Where n^, n/ are boson and fermion densities per square 
lattice site. The zero temperature fermion compressibil- 
ity, up to a Landau parameter correction, is approxi- 
mately equal to the Fermi pockets density of states, by 



(44) 



The boson compressibility, (using the xy model represen- 
tation of hard core bosons) is approximately 



32J fc 



(45) 



At zero temperature, ignoring boson-fermion interac- 
tions, we use Eqs.(|43|) and (|i^ ) to obtain the change in 
chemical potential to linear order in doping x 



(2K b + «/) 



(46) 



In the underdoped regime, where x « 1, the energy 
distance between \i and the fermion saddlepoints k sp ss 
(■7T, 0) defines the pseudogap A PQ as measured in tunneling 
and photoemmission (see Fig. fill. Its doping dependence 
is simply connected to the chemical potential shift 



A pg (x) = E k s P - n(x, T), 



(47) 



which yields a steady reduction of the pseudogap as a 
function of doping as plotted in Fig. [l2] 
In the normal state above T c , we have a theory of two 
decoupled, non interacting gases. For the bosons, we use 
a constant density of states and for the fermions 

we choose p/(e) from the dispersion (j36"l), (see Fig. [ll]). In 
this simplified theory, fi(T, x) can be found using Eq.(|4^) 
and solving 



e t — l 



dc- 



Pf( e 



e t +1 

The grand potential and entropy are given by 
fi(T) = -T [ dwp b {uj) log f 1 - e " ( ""'" (r> 



(48) 



+ T j dep f (e) log(e^ 
S(T) = -dn/dT. 
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FIG. 12. Thermodynamics of the PBFM. Excess nor- 
mal state hole entropy as a function of temperature for dif- 
ferent doping levels, calculated using Eqs. @. Inset: The 
pseudogap energy as a function of doping. Energies and tem- 
peratures are in units of the holes hopping parameter t' of Eq. 
©■ 



(49) 



In Fig. 12 the excess hole entropy S(T,x) of in the 
non superconducting state is shown for the density of 
states given in Fig|ll]. The picture which emerges is 
that above the superconducting transition temperature, 
bosons evaporate into the fermions gas. The evaporation 
is driven by the larger density of states of hole fermions 
than the bosons. This evaporation also implies a rapid 
increase in magnetic susceptibility. Its effects on trans- 
port have not yet been calculated. 



V. SUMMARY AND DISCUSSION 

This paper is primarily aimed at demonstrating the ap- 
plication of CORE to the Hubbard model, which allows 
us to extract its low energy degrees of freedom and derive 
the Plaquette Boson-Fermion model. The CORE calcu- 
lation could be improved by diagonalizing larger clusters 
within contemporary computational capabilities. A con- 
sistent computation of both the boson and fermion pa- 
rameters would be useful. It would permit systematic 
studies of extended Hubbard models and the effects of 
additional interactions. 

The d-wave hole pairs are already present in the Hubbard 
model on a single plaquette. Fortunately, due to the 
short coherence length and large hopping rate, the pairs 
maintain their integrity in the square lattice. 
The PBFM, at the simplest level of approximation, pro- 
vides a phase diagram which shares the basic features 
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of underdoped cuprates: the antiferromagnetic Mott in- 
sulator and a d-wave superconductor with nodal hole 
fermions. In the superconducting phase, the local spin 
one magnons are gapped at the antiferromagnetic reso- 
nance energy, and the remaining gapless excitations con- 
sist of a small density of hole pair bosons and spin half 
hole fermions. 

The PBFM brings us closer to understanding low tem- 
perature correlations of cuprates. It is amenable to mean 
field, low density, and variational approximations which 
do not lend themselves directly to the higher energy Hub- 
bard model and its various extensions. 
Here, the PBFM was only preliminarily explored. It 
would be interesting to study its thermodynamics and 
transport properties in more detail. 
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APPENDIX A: THE COMPLETE FOUR BOSON 
MODEL 

Here we present the complete four boson model including 
all interactions generated by CORE up to 2 plaquette 
terms. Coupling parameters are listed for square lattice 
and ladder geometries. We then estimate the magnitude 
of the truncated three plaquette terms. 
The four boson model can be separated into a bilinear 
part and a quartic part in the bosonic operators: 

n Ab = n b [b} + n t [t} + n mt [b,t] (Ai) 

where the bosons obey local hard core constraints 

&l&i+£4*m<l. (A2) 



The kinetic (bilinear) terms as written in section IV A 
are 



n 1 



(e 6 -2/i)^6t& i -J 6 ^(&J &i+ H.c 

Ju 



a{ij) 



H.c. 



(A3) 



The higher order interaction terms are 

[^o(ti*,-)o(*i*i)o 



7C 



V b 



2 , nun hj 



+V 1 (t l t 3 )\(t l t J ) 1 + V 2 {t l t J )\{t l t J ) 2 } 

-Ju x)( 6 iMi,-*«i+ii.c.) 



bt 

(ij) 



+Vbt V (plbitUaj + b]b/ m t m ), (A4) 



where (i»tj)s creates two triplets on plaquettes i and j, 
which are coupled into total spin S. When Vq — 2V\ = 
— 2V2 the triplet interactions may be written using spin- 
1 operators as • Sj. Similarly, For J t — Ju, which 
is close to the value given by CORE, (see table |), the 
bilinear two site triplet terms may be simplified to Jt^i ■ 
n 3 , with n a = ^(4 +t a ). 

The full Hamiltonian ^l] may serve as a starting point 
for various approximations or numerical studies. Its pa- 
rameters were computed using CORE from the Hubbard 
model with U/t = 3,6,8,10. The parameters are listed 
in tabic fl. 





U = 3t 


U = 6t 


U = 8t 


U = Wt 


eo 


-6.613 


-8.332 


-9.865 


-11.549 




(-6.019) 


(-7.983) 


(-9.593) 


(-11.324) 


«t 


0.152 


0.183 


0.174 


0.162 




(0.192) 


(0.263) 


(0.253) 


(0.233) 




1.178 


2.081 


3.557 


5.183 




(0.440) 


(3.212) 


(4.835) 


(6.567) 


Jt 


0.615 


0.397 


0.309 


0.249 


Ju 


0.590 


0.379 


0.297 


0.242 


Vo 


-0.361 


-0.152 


-0.114 


-0.099 


Vi 


-0.203 


-0.117 


-0.095 


-0.082 


V2 


0.214 


0.099 


0.071 


0.055 


Jt 


0.413 


0.340 


0.311 


0.289 


Jbt 


-0.383 


-0.233 


-0.173 


-0.134 


V u 


-0.133 


-0.286 


-0.143 


-0.191 


v bb 


0.884 


1.061 


1.145 


1.213 



TABLE I. Parameters for the Four Boson model, in units 
of t on the square lattice and ladder. The parameters were 
computed from the Hubbard model using range-2 CORE. Val- 
ues for the ladder are given in parenthesis where they differ 
from the square lattice. 

Note that the on-site terms for the ladder geometry 
(given in parenthesis in table |) differ from the square 
lattice case due to contributions of 2 plaquette terms 
hij . For example let e° be the bare on site triplet energy 
from the single plaquette spectrum and Set the correc- 
tion due to the inter plaquette interaction as described 
in section 



III A . The renormalized on-site energy at site 
i is €j + ZiSe t where Zj is the coordination number of site 
i. The values of and Set may be extracted from the 
table. For example: 
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Se t 
o 



{ square 

0* 



ladder 

4<5e t . 



)/2 



(A5) 



Estimation of the truncation error. In Fig. we com- 
pare between the low energy spectrum of the exact and 
the truncated effective Hamiltonian for 3 collinear pla- 
quettes. This comparison may be used to estimate the 
magnitude of the higher order three plaquette terms hijk 
defined by equation (^Tj) 



h 



rjk 



Tjren 



h 



jk 



hi + hj 



rrren 
n ijk 



H eff 

(A6) 



Recall that H^J} has the exact low lying spectrum of the 
original Hamiltonian on the 3 plaquettes. Thus expecta- 
tion values of h^k in the ground state and first excited 
states are calculated by subtracting energies of Hff^ from 
corresponding exact energies of the three plaquette prob- 
lem. We estimate (hij) in a similar way, by comparing 
energies of two disjoint plaquettes to the exact energies of 
two coupled plaquettes. Small expectation values (hijk) 
relative to (hij) suggest rapid convergence of the cluster 
expansion. 





S=0 


S=l 


holes 


330 


7.7 


2 holes 


27 


19.5 



TABLE II. Convergence of the cluster expansion. 

The ratio (hij) / (hijk) given for different sectors in the Hamil- 
tonian with U = &t indicates excellent convergence of CORE 
on a ladder. 



APPENDIX B: CORE CALCULATION FOR THE 
TIGHT BINDING MODEL 



In section III B we discussed the effects of breaking lattice 
transnational symmetry within the reduced Hilbert space. 
We argued that interactions of increasing range gradually 
reduce the effects of symmetry breaking on the spectrum. 
It is instructive to study this process in a simple model 
where a CORE calculation can be carried easily to long 
ranges. Such an opportunity is provided by the tight 
binding model on a chain: 



x Ci). 



(Bl) 



E/t 



-22 
-22.2 
-22.4 
-22.6 
-22.8 

-23 
-23.2 
-23.4 
-23.6 
-23.8 

-24 



S=1 

s=o 



S=1 

s=o 



o 

CM 



c 



Effective H 



Exact 



FIG. 13. Low energy Spectrum of exact compared to 
effective Hamiltonian on 3 plaquettes. The comparison 
is presented for the Hubbard model with U=6t in the 0-hole 
and 2-hole sectors. An arbitrary chemical potential was used 
to set the 2-hole energies slightly above the plotted 0-hole 
energies. 

Table H gives a summary of the ratios (hij) / (hijk) in the 
lowest states of the different sectors of the Hamiltonian. 
The satisfactory convergence of the cluster expansion, 
implies the integrity of bosonic states on the lattice, at 
least for ladder geometry. Interestingly, it is still very 
good for the Hubbard model with U = 6t where pair 
binding energy is positive on a plaquette. This strength- 
ens the argument that binding is generated dynamically 
on the lattice. The holes remain tightly bound because 
correlated motion reduces their kinetic energy. 



We apply CORE to the single electron sector of this 
model, coarse graining it to blocks of 2 sites. In each 
block we retain only the empty state | ) i and the single 
electron symmetric state of energy — t: 



//|0) 



1 



V2 



(c 



4i+l) 



(B2) 



Hence we can only hope to reconstruct the lowest of the 
2 bands in the folded Brillouin zone k — [— tt/2, ir/2}. 
The effective Hamiltonian generated by CORE at any 
range of the cluster expansion is of the general form: 



(B3) 



Such a Hamiltonian cannot reproduce the sharp band 
edge at k = ±7r/2 at any finite range of hopping. How- 
ever, as demonstrated in figure (p^), CORE calculations 
of increasing range introduce higher harmonics that suc- 
cessively approximate the sharp edge. If one is interested 
in the properties of the model far from the dimerized zone 
edge then by Fig, 14 the effective Hamiltonian generated 
by range-3 CORE should suffice. 
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